Failure of the Wiedemann— Franz Law in Mesoscopic Conductors 
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We study the effect of mesoscopic fluctuations on the validity of the Wiedemann-Franz (WF) law 
for quasi-one-dimensional metal wires and open quantum dots. At temperatures much less than 
the generalized Thouless energy, E c , the WF law is satisfied for each specific sample but as the 
temperature is raised through E c a sample-specific correction to the WF law of order 1/g appears 
(g is the dimensionless conductance) and then tends to zero again at ksT S> E c . The mesoscopic 
violation of the Weidemann- Franz law is even more pronounced in a ring geometry, for which Lorenz 
number exhibits h/e flux-periodic Aharonov-Bohm oscillations. 

PACS numbers: 73.23.-b, 72.15.Eb, 73.63.Kv 
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I. INTRODUCTION 

The Wiedemann-Franz (WF) lawi in macroscopic con- 
ductors relates the electrical and thermal conductivities 
a, k via the condition n/aT = l , where the Lorenz 
number lo = (n 2 /S)(kp,/e) 2 (here k-g,e are Boltzmann's 
constant and the electron charge). The law is satis- 
fled with an accuracy of 5 — 10% in most metals at 
room temperature 2 , but fails in the temperature range 
T w 10 — 100K where the Bloch-Gruneisen form of the 
resistivity holds^. In this range the failure is under- 
stood because inelastic phonon scattering degrades ther- 
mal and electrical currents differently. At temperatures 
below 10K, when the inelastic contribution to the scatter- 
ing is negligible the law is satisfied extremely accurately^ 
as it follows simply from the Sommerfeld expansion of 
the energy-dependent conductivity, which is accurate to 
order (k-BT/Ep) 2 , where Ep is the Fermi energy. Fail- 
ure of the WF law at low temperatures is often asso- 
ciated with non-Fermi-liquid behavior in strongly corre- 
lated systems^. 

In the current work we focus on the violation of the WF 
law induced by mesoscopic fluctuations in Fermi-liquid 
systems. Mesoscopic conductors are normal conductors 
on the micron or nanometer scale connecting bulk leads 
which act as thermal equilibrium reservoirs held at differ- 
ent voltages and temperatures. By assumption there are 
no inelastic scattering processes in the conductor itself, 
exactly the condition under which the WF law holds for 
macroscopic conductors. However, as is well known, each 
specific mesoscopic sample exhibits a fluctuating energy- 
dependence of its conductance which varies by SG ~ e 2 /h 
on an energy scale E c <C Ep&l. Therefore the Sommer- 
feld expansion fails even when k^T <C Ep. The most 
immediate effect of this failure is that the thermopower 
coefficient, which gives the voltage induced by a tem- 
perature gradient when no electrical current flows, be- 
comes random in sign and is enhanced by (Ep / E^g^ 1 , 
where g = G/Gq is the dimensionless conductance and 
Go = e 2 //i. This enhancement can be many orders of 
magnitude, has been discussed in a number of previous 
works 8,9,10,11 - 12 , and has been measured in at least one 
experiment on quantum dotsi^. 



It is straightforward to show that the violation of the 
WF law is not due to the enhanced thermopower coeffi- 
cient and occurs at size scales for which the thermopower 
effect would be negligible. As usual, the mesoscopic cor- 
rection to the WF law is random and arises from inter- 
ference of multiply-scattered electron trajectories; in a 
ring geometry this correction will oscillate periodically 
in magnetic flux with period h/e, making it easier to 
detect. Very recently thermoelectric measurements have 
been made in such systems^, so experiments of this type 
are feasible. 

We consider a two-probe measurement for which the 
electrical current, I c , and thermal current, 1^, between 
the reservoirs at temperatures T and T + 8T and at volt- 
age V are 



I e = GV + BST, 
I h = TV + EST. 



(la) 
(lb) 



The off-diagonal coefficients, B and V, are related by 
Onsager's principle: T — —TB. 

To the leading order in (Epr) -1 the electron-electron 
interaction does not affect the fluctuations of transport 
coefficients at low temperature (r is the mean elastic 
scattering time)i&; in addition inelastic scattering due 
to interactions at finite temperature may be neglected in 
the mesoscopic regime. Interactions in disordered metals 
do affect the average electric and thermal conductances 
and can lead to a violation of the WF law^i^; however 
this effect does not have magnetic flux sensitivity and 
its temperature dependence is different from that of the 
mesoscopic violation studied here. 

Neglecting the interaction effects, the thermoelectric 
coefficients can all be expressed in terms of integrals 
over G(e), where the T = conductance at energy 
e is given by the Landauer-Buttiker formula G(e) = 
2Go Yl a 8 (£■) | 2 and t a p(e) is the transmission ampli- 
tude for conduction channels a, (3. Defining the energy 
moments of the conductance G(e) at temperature T and 
chemical potential \x by 



G {n) {T) 



; n (#,4(e), G(£ )M| ( 2) 
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FIG. 1: Upper panel shows the thermal weight functions for 
the electrical current, I c , and for the heat current, Ih. For 
constant electron transmissivity ~ G(e) on the scale of ksT 
the ratio of areas below these functions leads to the WF law, 
which thus always holds as T — > 0. However if SG(e) fluc- 
tuates on the scale of E c ~ ksT (lower panel) then 8G(e) is 
unequally weighted for I e and Ih and sample-specific correc- 
tions to the WF law appear. 



one finds that G = G^°\ B = -{k B /e)G^ and 5 = 
{k B /e) 2 TG {2 \ The thermal conductance G h = h/ST, 
measured when 7 e = 0. From Eq. Q we have Gh = 
E — TB/G. The WF law has the same form for these 
conductance coefficients as for the bulk conductivities: 
Gh/GT — l . 

If G(e) is constant on scales much greater than k B T 
then the thermal conductance Gh obeys the WF law to 
very high accuracy. First, if G(e) is constant and hence 
equal for electrons and holes within k B T of the Fermi en- 
ergy then no net electrical current flows between the two 
reservoirs under a thermal gradient and the thermopower 
is zero: T = —TB = 0. Thermal current does flow since 
the colder reservoir receives more electrons with energy 
e — /i > and holes with energy e — fi < 0; this flux is 
given by setting G(e) = G((i) = G M , leading immediately 
to G& = (tt 2 /3)G m and G h /G = E/G = l T. Thus the 
WF law does hold as T — > for mesoscopic conductors; 
fluctuations in both thermal and electrical conductances 
maintain a fixed ratio given by the WF law. 

If G(e) does vary on the scale of fc^T then one can- 
not approximate G(e) by G M , the thermal and electrical 
currents have different weighting factors (see Fig. 1), 
and their ratio is no longer fixed at the WF value. Thus 
mesoscopic conductors will begin to violate the WF law 
at temperatures which are on the scale of variation of 
G(e). This scale is E c — h/r t , where r t is the typical 
transit time between the reservoirs and depends on the 
geometry of the conductosSi Since the mesoscopic con- 
ductance fluctuates by 8G ~ e 2 /h on this scale we can 
immediately estimate the maximum size of the violation 
of the WF law. Define the Lorenz ratio to be 
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5 TB 
GT ~ TCP'' 



(3) 



fluctuating contributions SG(T), 5E(T) are ~ g l ; there- 
fore the fluctuating part of the first term in Eq. © is 
5(E/G) = l Q [8E/l Q G^T - <5G/G P ] - g-\ In contrast, 
the thermopower coefficients T, B are zero in leading or- 
der and only have a contribution from the mesoscopic 
fluctuations2*2iii; therefore the fluctuating part of the 
second term in Eq. © is TB/G 2 ~ g~ 2 and the ther- 
mopower fluctuations do not contribute to the violation 
of the WF law at leading order. Thus the fluctuations 
lor] in the Lorenz ratio I are given by 



l = h{l + ri) = 



GT' 



1 

G„ 



SE 
l a T 



-5G 



(4) 



The standard deviation of the dimensionlcss quantity 
r] gives us the typical violation of the WF law. Be- 
fore averaging for k B T <C E c we can expand G(e) in 
Eq. (J2J around // = E-p to order (e — fi) 2 ; the zeroth 
order yields the WF law, the first order vanishes by 
symmetry, and the second order will give a correction 
7] ~ {k B T) 2 G"{^)/G^, which for a mesoscopic sample 
will be <~ (k B T / E c ) 2 g -1 . Hence the deviation from the 
WF law will rise quadratically for fcgT -C E c and become 
— g- 1 at k B T — E c . When k B T > E c the typical vi- 
olation will again decrease as the random contributions 
from energies e begin to self-average over the window 
kgT. The quantitative behavior of the fluctuations can 
be obtained from the variance of r), which is determined 
by the zero temperature conductance correlation func- 
tions K(( £ - £')/E c ) = (SG{£)SG{£ r ))/G 2 as 



var(?7) 



r* 2 

(G,) 2 



9G 2 



6G 2 



+ G, 



on 



f 2 (T,E c ) 
9 2 



(5) 

where the coefficients C nm = (SG^SG^) jG\ refer to 
integrals of the form: 



G nm — / d£al£ 



jdj df (£ - fl) n (£' - HY 

d£ cfe' (k B T) n 



+m 



■K 



E c . 



(6) 



where G = G p + SG(T) and 3 = ZqTG^ + 6E(T). The 



As already noted, the typical value of the violation of the 
WF law is rms(?7) = f(T,E c )/g; the function f(T,E c ) 
is plotted for the geometry of a quasi-one-dimensional 
metallic wire and an open quantum dot in Fig. 2. As 
expected, the maximum of f(T,E c ) occurs at k B T ~ 
E c and /(T, E c ) vanishes for both k B T <§C E c and for 
k B T ^> E c . We now provide a more detailed analysis of 
the results. 



II. THERMAL CONDUCTANCE 

A. Quasi-one-dimensional diffusive conductors 

For a diffusive mesoscopic system (i.e., a system 
with the elastic mean free path shorter than all sam- 
ple dimensions) the conductance correlation function 
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FIG. 2: Temperature dependence of the rms(ry) in quasi-one- 
dimensional wires (solid line) and in open quantum dots with 
N\ = iV r 3> 1 (dashed line). In both cases a broad maximum 
in the violation of the WF law occurs at fcsT ~ E c and 
persists to much higher temperatures. 



K(Ae/E c ) has been calculated by impurity-averaged per- 
turbation theory§i£ and can be expressed in terms of 
sums over the eigenvalues of the diffuson and cooperon 



propagator; 



,18.19 



For the case of a disordered quasi-one- 



dimensional wire of length L with diffusion coefficient D, 
we have E c = D j L 2 and 



K{x) 



71—1 



2x 2 



e+x 2 (e+x 2 ) 2 ' 
. (7) 

where the parameter /3 = 1 in the presence of time- 
reversal symmetry and (3 = 2 in its absence, and we 
have included spin degeneracy. Time-reversal symmetry 
is effectively broken by a magnetic field comparable to 
a flux $o = through the wire, typically rather weak 
fields. Experiments which rely on varying magnetic field 
to observe fluctuations will correspond to the case (3 = 2. 
The function f(T, E c ) giving rms(ij) for a wire is plotted 
as the solid line in Fig. 2; it has the asymptotes: 



E, 



f{T,E c 



A_ An 2 - 30 [ET 
3^V 5^ V fenT 



; k B T < E c , (8a) 
/c B T> B c .(8b) 



The maximum of f(T,E c ) is f max {T,E c ) » 0.23 at 
^BTmax ~ 0.5E C leading to a violation of the WF law 
~ 0.23/(7- This expression describes both micron scale 
and nano scale metallic wires as well as semiconducting 
quantum wires in the diffusive regime for j » 1. For 
such systems the dimensionless conductance g can vary 
between 10° — 10 3 . When g approaches unity it will be 
necessary to include the thermopower corrections to the 
WF law as well and the violation, while order unity, will 
not be quantitatively described by f(T, E c ). 



B. Aharonov— Bohm oscillations 

In order to make the mesoscopic deviation from the 
WF law more easily measurable it will be convenient to 



fabricate ring structures similar to those used to first 
measure the mesoscopic /i/e-periodic Aharonov-Bohm 
(AB) effect in normal metals 20-21 . In this case the meso- 
scopic quantum corrections to the conductance are the 
only terms which give rise to oscillations of the conduc- 
tance with period h/e and these terms are relatively 
easy to isolate by Fourier transforming the measured 
behavior of G(<&) and Gh($) ~ 2($), where $ a 5 
is the magnetic flux through the ring in magnetic field 
B. We consider then the ratio of harmonics of the AB 
oscillations at frequency k/$o,k = 1,2,... in <5G( ( &) 
and <5S(<&). For k-gT -C E c we have the simple re- 
sult that <SG($) ra <5G M (*) and 53 » Z T<5G M ($); hence 
<5S($)/£G($) = l Q T for all $. The AB oscillations are 
perfectly in phase and their ratio is exactly the Lorenz 
number (as is the ratio of each harmonic). However 
for fegT ~ E c the oscillating contributions to G, H are 
no longer proportional to one another and their ratio 
need not be the Lorenz number. To quantify their rela- 
tionship we calculate the harmonics of the conductance 
and heat conductance correlation functions (in flux) for 
such a ring. Let F G (A$) = (<5G($)<5G($ + A$)) and 
Fh(A$) = (<5S($)(5E:($ + A$)); neglecting the flux 
through the wires making up the ring (i.e. including only 
the flux through the hole) these functions are periodic 
with fundamental frequency $0 an d can be represented 
in terms of their harmonics F, 



(fc) 



F G , H (A$) = a^ fc=o 4 fe Uos 



/2tt£;A$\ 



(9) 



where a is geometry dependent factor. As above, Eq. J7J), 
the relevant correlation functions can be expressed as 
sums where for the ring the integers in the sum appear 
only in the combination m — A$ / $0 (leading to an ob- 
vious periodicity in <&o). It is convenient in this case 
to rewrite the sum using the Poisson summation for- 
mula where the integration variable is shifted to <f> = 
to — A$/$q 21 . This yields a relatively simple formula for 



for the harmonics F G = GqF 00 of the electrical con- 



ductance and for the harmonics F^ = (k^T / h) 2 F^' of 
the thermal conductance. Here 



de de' 



(k B T) 



d(/>cos(27rfc0)/C U 2 + -f-, 



7 e — e 



Er E r 



(10) 



with kernel /C(£, x) defined by Eq. (7J|, Er = D/R 2 is in- 
versely proportional to the electron diffusion time along a 
circumference of the ring of radius R, and 7 takes into ac- 
count electron escape from the ring through the leads and 
other phase relaxation rates, such as electron-electron or 
electron-phonon dephasing rates. 

Performing the integrations in Eq. (|10fl . we find that 
the ratio of the harmonics of S to G decreases from the 
WF value monotonically with increasing temperatures 
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FIG. 3: (Top) Temperature dependence of Fg~ > (T) for two 
values of "//Eh- At high temperatures ksT ^> Er amplitudes 
f g,b 1] ( T ) decrease only as T" 1 . (Bottom) the ratio /f£ ] 
as a function of temperature. At T > Er the ratio saturates 
to 0.55(foT) 2 (dotted line), see Eq. ATI . 



until it saturates at k B T > En, 7, to the value 



1 F 



(fc) 



T 2 p( fc ) 



21n 2 - 180 

57T 2 



li « 0.55^ 



(11) 



In this temperature regime the conductance and heat 
conductance AB oscillations will appear uncorrelated and 
their ratio will randomly vary at each value of magnetic 
flux $; it is only the ratio of their variances (or harmon- 
ics) averaged over magnetic fields which will saturate to a 
constant value. This saturation value will be different by 
roughly a factor of one half from the value expected from 
the WF law. The reason that the saturation value is less 
than the WF value when k B T > E c is that the conduc- 
tance fluctuations are determined by transmission fluctu- 
ations within k B T of the Fermi level, while the thermal 
conductance fluctuations are determined by transmission 
fluctuations in two regions of width k B T displaced from 
the Fermi level by ~ ±k B T, see Fig. 1. Hence the con- 
tributions from different energies to the thermal conduc- 
tance fluctuations are less correlated and tend to cancel 
each other more than the contributions to the electrical 
conductance fluctuations. We plot the ratio /F(p as 
a function k B T '/Er for several values of 7 in Fig. [31 



C. Open Quantum Dots 

If we consider semiconducting systems the above cal- 
culations also apply and since g is typically much smaller 
one expects large violations of the WF law, although in 
general the WF law is not well-satisfied in semiconduc- 
tors so the novelty is less than in metals. Semiconduct- 
ing quantum dots exhibit large aperiodic thermopower 
fluctuations as a function of magnetic fielcU^; these have 
been treated theoretically by Van Langen et aliii. The 
behavior of the Lorenz number for this case was not stud- 
ied either experimentally or theoretically. As long as the 



FIG. 4: Temperature dependence of the rms(B) in quasi-one- 
dimensional wires (solid line) and in open quantum dots with 
N\ = N t 1 (dashed line). The maximum in rms(B) occurs 
at fc_gT ~ E c and vanishes at low (fcsT -C E c ) and high 
(fcflT 2> E c ) temperatures. 



number of open channels W>1 our Eq. (JjJJ for the vari- 
ance of rj still applies with the only change being the 
replacement of Eq. for the conductance energy corre- 
lation function by the appropriate function for an open 
quantum dot. If we assume the quantum dot junction 
generates chaotic scattering then the S-matrix of the sys- 
tem will be described approximately by the Dyson ensem- 
bles of random matrix theory and the energy correlation 
function K qd (Ae/E c ) = (N?N? /N* h )(l + Ae 2 /E 2 )~\ 
where E c — N c \ l Si/(2n), N\, N r are the number of incom- 
ing channels in the left and right leads of the junction, 
A'ch = N T + N\, and 61 is the mean level spacing of elec- 
tron states in the do1s22. Substituting K qc \(x) into Eq. JJjJ, 
we calculate var(?y) = f* d (T, E c )/g 2 with g = 2NiN T /N ch 



from Eq. (JSJ 
forms: 



f q d(T,E c ) has the following asymptotic 



/qd 
/qd 



16tt 2 2N,K 



3 iVi 



k B T 



E r 



2tt 2 - 15 N\N T 



15tt 



E c 
k B T 



k B T < E c , (12a) 
k B T » E c , (12b) 



and reaches its maximum value /^ ax ~ 0.88NiN T /N^ h at 
k B T max w 0.68.E C . We plot the function f qd (T,E c ) for a 
chaotic quantum dot (dashed line) in comparison to the 
wire case in Fig. 2; the behavior of the two cases is quite 
similar. Note that in both cases the decay at k B T ^> E c is 
very slow ~ J 1-1 / 2 and the maximum is very broad on the 
high temperature side. Finally we note that for a non- 
chaotic quantum dot or metallic nanobridge 23 the con- 
ductance fluctuations can exceed Go and consequently 
the violations of the WF law can be even larger than we 
find here. The behavior for specific non-chaotic shapes 
is non-universal and can be calculated by semi-classical 
methods^. 



III. THERMOPOWER 

For completeness we make a few brief comments about 
the thermopower coefficient, which has been treated in 
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RefsAS for the case of a disordered wire and in Refiii for 
the case of a chaotic quantum dot. The thermopower 
Q = —B/G is the coefficient of the voltage induced 
by a temperature difference when no electrical current 
flows, I c = 0. To leading order in the Sommcrfcld ex- 
pansion Q = and at first order Q ~ {k B T '/e)G^/G M ~ 
(k B T / E-p) for a macroscopic conductor. For a mesoscopic 
conductor G^/G M ~ (l/gE c ) and the thermopower fluc- 
tuates in sign and for k B T ~ E c is enhanced by a factor 
{Ey / E^g -1 which is typically several orders of magni- 
tude. The fluctuations of Q to the lowest order in 1/g 
are determined by fluctuations of B, the latter can be 
calculated from Eq. JBJ) for C\\\ 

var(B) = ft^l\ 2 0u ee (*°£±\ 2 h 2 (T, E c ) . (13) 

The behavior of rms(B) = (k^G^/ e)h(T, E c ) as a func- 
tion of k-sT/E c is similar but not identical to that of 
var(77), increasing linearly with T for k B T <C E c and then 
decaying as T -1 / 2 for k B T » i? c , see Fig.0] For diffusive 
wires /i(T, E c ) was calculated in RefiS. Here we present 
the low and high temperature asymptotes of h(T,E c ): 

2 k T 

h{T,E c ) = -VlOC(8)-|-, fc B T«S c , (14a) 

fc(T,*0 = V^VS' fcBT>Sc (Wb) 

with £ c = D/L 2 . The maximum of h(T,E c ) is /i max rs 
0.17 at ~ 0.37i? c . For geometry of an open quantum 
dot, E c = N c hSi/(2w) and the asymptotes are 

h qd {T,E c ) = — -2 — , k B T«.E Cl (15a) 

6 M ch E c 

ftqd( r, £c )=^E^^^, te r» £o . (1 5b, 



The maximum of h w OSGNiNz/N^ occurs at fceT ~ 
0.49i? c . The low temperature asymptote, Eq. I|15a|) . was 
previously found in Ref . IllL 

The maximum fluctuations in thermopower are 
Qmax ~ {^B/s)g~ 1 at k B T oj E c . The thermopower 
corrections to the thermal conductance Oh are always 
smaller than the direct quantum correction by a factor 
g _1 and therefore only play a significant role in the vio- 
lation of the WF law when g approaches unity. 



IV. CONCLUSIONS 



In conclusion, we evaluated the variance of the ratio 
of the heat and electric conductances in mesoscopic sys- 
tems. We showed that this ratio does not fluctuate as 
T — > and exactly satisfies the Wiedemann-Franz law, 
but fluctuates with a finite variance for non-zero T and 
is a non-monotonic function of temperature with a max- 
imum at temperature equal to the Thouless energy. The 
ratio of the e/h harmonic of the heat and electric con- 
ductance correlation functions in an Aharonov-Bohm ge- 
ometry satisfies the Weidemann-Franz law only at low 
temperatures; at high temperatures their ratio becomes 
0.55Z 2 . 
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M. Keck Foundation and by NSF Materials Theory grant 
DMR-0408638. 
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